The crystal structure of cold compressed graphite 
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Through a systematic structural search we found an allotrope of carbon with Cmmm symmetry 
which we predict to be more stable than graphite for pressures above 10 GPa. This material, which 
we refer to as Z-carbon, is formed by pure sp^ bonds and is the only carbon allotrope which provides 
an excellent match to unexplained features in experimental X-ray diffraction and Raman spectra of 
graphite under pressure. The transition from graphite to Z-carbon can occur through simple sliding 
and buckling of graphene sheets. Our calculations predict that Z-carbon is a transparent wide band 
gap semiconductor with a hardness comparable to diamond. 
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Thanks to the flexibility to form sp, sp^ and sp'^ bonds, 
carbon is one of the most versatile chemical elements. 
At ambient pressure, it is usually found as graphite (the 
most stable structure) or as diamond, but the richness of 
its phase diagram does not end there. In fact, many other 
structures have been proposed during the past years, es- 
pecially since experimental data suggested the existence 
of a super hard phase of carbon. Evidences for a struc- 
tural phase transition in compressed graphite to this un- 
known phase of carbon have been reported in numerous 
experiments [H-Q]- fact, in the range of 10 to 25 GPa 
one observes an increase of the resistivity [J and of the 
optical transmittance 0, Q : a marked decrease of the o^ 
tical reflectivity changes in near fc-edge spectra 
and in X-ray diffraction (XRD) patterns [5|-l7[. Several 
hypothetical structures have been proposed to explain 
these features, such as hybrid sp^-sp^ diamond-graphite 
structures [1], M-carbon [§], bct-C4-carbon [l3| and W- 
carbon However, none of these structures is able to 
match all experimental data in an unambiguous and fully 
satisfactory manner. 

A common way to search for new crystal structures is 
to perform a systematic survey of the enthalpy surface 
using some sophisticated structure prediction method 
(for discussion on such methods see Ref. [13 )• Here we 
use the minima hopping method [l^ (MHM) for crys- 
tal structure prediction [l^ . which was designed to ex- 
plore low-enthalpy phases of materials. This method 
was coupled to the all-electron projector-augmented wave 
method as implemented in the ABINIT code [l^, [l6| . 
Within the MHM, the system is moved from one config- 
uration to the next by performing consecutive molecular 
dynamics escape steps and geometry relaxations. The 
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FIG. 1. (a) Structure of Z-carbon viewed from two differ- 
ent angles revealing planar four-membered and non-planar 
eight-membered rings forming chains along the fe-direction 
and channels in the c-direction. The graphene sheets are in 
the a-c plane, (b) Proposed transition pathway from graphite 
to Z-carbon. 



initial velocities for the dynamics are aligned preferably 
along soft-mode directions in order to favor the escape 
to low enthalpy structures. Revisiting already known 
structures is avoided by a feedback mechanism. Relax- 
ations are performed by the fast inertia relaxation en- 
gine 



17| . The local density approximation was employed 



based on its good description of graphite. However, the 
enthalpy ordering was reconfirmed within the general- 
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izcd gradient approximation using two different function- 
al (PBE Q and PBEsol (T^l). The most promising 
candidate structures were tfien re-relaxed using norm 
conserving Hartwigsen-Goedecker-Hutter pseudopoten- 
tials [l^l- Carefully converged Mankhorst-Pack fc-point 
meshes were used together with a plane wave cut-off en- 
ergy of 2100 eV. 

The MHM was employed using simulation cells with 
4 and 8 carbon atoms at a constant pressure of 15GPa. 
We found, in addition to previously proposed structures 
of cold compressed graphite, a carbon phase that we 
call Z-carbon. This structure has Cmmm symmetry (see 
Fig. [1^) and, like diamond, is composed of sp^ bonds. 
The conventional unit cell has 16 atoms with cell pa- 
rameters at OGPa of a = 8.668 A, h = 4.207 A, and 
c = 2.486 A, yielding a cell volume of Vb = 90.7 A^. 
The two inequivalent carbon atoms occupy the 8p and 
8(7 crystallographic sites with coordinates (1/3, y,0) and 
(0.089,2/, 1/2), where y = 0.315. The structure contains 
four-, six- and eight-membered rings, where planar four- 
membered rings and non-planar eight-membered rings 
join together buckled graphene sheets. This structure 
can be interpreted as a combination of hexagonal dia- 
mond and bct-C4-carbon plj . 

In contrast to other structure prediction methods like 
evolutionary algorithms or random search, the efficient 
escape moves in the MHM are based on fundamen- 
tal physical processes. Therefore, minima found con- 
secutively during a MHM simulation are usually con- 
nected through low enthalpy barriers. Since we have 
observed escape moves to and from Z-carbon to occur 
exclusively from and to graphite, we expect this tran- 
sition to be the most probable. In Fig. [T}) we show a 
possible transition pathway from graphite to Z-carbon. 
This process is a combination of sliding and buckling of 
the graphene sheets. The naturally staggered, i.e. AB 
stacked, graphene sheets slide along the [210] direction to 
an aligned AA stacking while the inter-layer distance de- 
creases, and the aligned graphene sheets deform to create 
an alternating armchair-zigzag buckling. 

In order to investigate the relative stability of Z- 
carbon, the calculated enthalpy difference with respect 
to graphite of several allotropes are compared in Fig. [2] as 
a function of pressure. Z-carbon has the lowest enthalpy 
among all proposed cold-compressed graphite phases, be- 
coming more stable than graphite at 9.9 GPa (around 
2.5 GPa below M^-carbon). 

We further investigated the dynamical lattice stabil- 
ity of this phase by computing the phonon dispersion in 
the whole Brillouin zone. We used linear-response the- 
ory in the framework of density functional perturbation 
theory [22j with the ABINIT code. A proper convergence 
was ensured with a 12x12x12 fc-point sampling, a 4x4x4 
g-point mesh, and a cut-off energy of 800 eV. All phonon 
modes were real confirming the strutural stability of the 
structure. Furthermore, from a fit of the Murnaghan 
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FIG. 2. Calculated enthalpy dilference per atom with re- 
spect to graphite of several carbon allotropes as a function of 
pressure. Graphite is the horizontal line at zero. Z-carbon 
becomes more stable than graphite at around 10 GPa. 



equation we obtained a bulk modulus of = 441.5 GPa, 
and using the method proposed by Gao et al. (23| we 
calculated a Vicker's hardness of H^, = 95.4 GPa. Both 
bulk modulus and hardness are extremely high and very 
close to the values for diamond (B^iamond ^ 463.0 GPa 
and ijdiaraond ^ 97.8 GPa), which is compatible with the 
observed ring cracks in diamond anvil cells . 

To investigate the energy gap of this material we used 
the perturbative many-body GW technique starting from 
the local density approximation These calculations 
reveal that Z-carbon is an indirect band-gap material 
with a gap of around 4.7 eV. Therefore, this material is 
expected to be optically transparent in agreement with 
experiments [3, 

We have gathered a considerable amount of experimen- 
tal evidence that points to the presence of Z-carbon in 
cold compressed graphite samjjles. The first comes from 
the XRD experiment of Rcf. [71 . In Fig. |3] we can see 
that the broadening of the XRD-spectra at high pres- 
sure can be explained by the coexistence of graphite and 
Z-carbon. However, the experimental curve can also be 
explained to some extent by the other proposed carbon 
allotropes so that this experiment alone is not con- 

clusive. 

Other strong signatures for Z-carbon can be gathered 
from our measurements of Raman spectroscopy under 
pressure. These experiments were carried out at 300 K 
using the 514.5 nm line excitation of an Ar"*" laser, and a 
Jobin-Yvon HR-800 Labram spectrometer with double- 
notch filtering with resolution better than 2cm^^. In 
the high pressure Raman measurements, we used a dia- 
mond anvil cell to apply pressure on two different sam- 
ples (single crystals of graphite and highly oriented py- 
rolitic graphite), inside a 120 micron hole drilled in an 
iconcl gasket. Argon and paraffin was used as the pres- 
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FIG. 3. Experimental XRD for cold compressed graphite at 
two difTerent pressures from Ref. 0] and simulated XRD pat- 
tern for Z-Carbon (at 23.9 GPa) and Graphite (at 13.7 GPa). 
The main characteristics of the proposed carbon are perfectly 
in agreement with the experimentally observed changes. 



sure medium. The pressure was determined by the ruby 
luminescence of a smah chip (< 30 microns). The laser 
was focused down to 3 microns with a power of about 
20 mW on the sample. 

The principal Raman active mode of graphite is the G- 
band at 1579 cm~^ (at OGPa) which originates from the 
sp^ carbon atoms vibrating in-plane with E2g symmetry. 
The effect of hydrostatic pressure on the linewidth of the 
G-band is shown in Fig. U] The linewidth remains nearly 
constant until around 9-10 GPa. Above this value, the 
linewidth begins to broaden rapidly, in agreement with 
previous results of Hanfland et al. [3| . (A similar broad- 
ening has also been reported for turbostratic graphite- 
like BC4 under pressure [1^.) This behavior is a sign of 
a structural transformation at this pressure, and can be 
explained by important changes in the Raman cross sec- 
tion caused by interlayer coupling and the formation of 
sp^ bonds. As seen in Fig.[2l Z-carbon becomes enthalpi- 
cally favored with respect to graphite at around 10 GPa, 
whereas all other proposed structures cross the graphite 
line at significantly higher pressures. 

Finally, we present the strongest direct evidence of 
the existence of Z-carbon, which is found in the Raman 
spectrum of graphite under hydrostatic pressure, shown 
in Fig. [5] for the energy range below the 1** order Ra- 
man peak of diamond (1332 cm^ at OGPa) Neither 
graphite nor cubic-diamond have Raman active peaks in 
the selected energy region, however we can observe that 
a clear peak appears at 1082 cm^^ for pressures higher 
than 9.8 GPa. This peak can not be explained by either 
bct-C4 carbon, M-carbon, nor by the pressure medium 



FIG. 4. Experimental linewidth of the G-band of graphite un- 
der pressure. The linewidth stays nearly constant until pres- 
sures of the order of 9 — 10 GPa, above which the linewidth be- 
gins to broaden rapidly. This is a strong evidence for a struc- 
tural transition in graphite. Experiments were conducted us- 
ing highly oriented pyrophillitic graphite (HOPG) and argon 
(squares) or paraffin oil (triangles) as pressure transmitting 
media. The black dots are taken from Ref. 01 



(argon). The only structures that have Raman active 
modes compatible with this experimental evidence are Z- 
carbon and W-carbon. For Z-carbon the frequencies are 
1096.5 cm-i at 10 GPa and 1110 cm"! at 15 GPa. In- 
cidentally, Z-carbon also has a Raman active Ag mode 
at 1348.5 cm^^ at OGPa (theoretical value). This ap- 
pears as a signature of planar four-membered rings that 
overlaps with the so-called defect D-band of graphite at 
around 1345.5 cm~^ at OGPa (experimental value). 

In conclusion, we identified an allotropic structure 
of carbon, Z-carbon, that becomes more stable than 
graphite above 10 GPa. From all known carbon al- 
lotropes, only cubic and hexagonal diamond have lower 
enthalpy at high pressures. The Z-carbon structure is as 
hard as diamond, and is transparent in the optical re- 
gion. Moreover, a wide range of experimental data can 
only be explained by the presence of Z-carbon in samples 
of cold compressed graphite: first, the features of the X- 
ray diffraction spectra of graphite under pressure exhibit 
a broadening that matches the main peaks of Z-carbon. 
Second, the principal Raman signal of graphite, the G- 
band mode, suffers an abrupt increase of the linewidth 
above 9-10 GPa — the pressure range where Z-carbon be- 
comes more stable than graphite. Third, a new peak at 
1082 cm^^ appears in the Raman spectrum of graphite 
at around 10 GPa, at the frequency of a Raman active 
mode of Z-carbon. The only candidate among all car- 
bon allotropes proposed so far that can explain all above 
features simultaneously is Z-carbon. 
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FIG. 5. Experimental Raman spectra of graphite under pres- 
sure. The peak around 1082 cm~^ appearing at around 
10 GPa and its evolution can be explained by either Z-carbon 
or W-carbon. 



Our work also highlights the promising prospects of 
the minima hopping method for crystal structure predic- 
The exploration of the structural variety of 
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tion 

even simple elements such as carbon was up to now typ- 
ically the subject of many different studies which were 
presented in numerous papers over many years. In this 
first application of the MHM we were able to find not 
only Z-carbon, but also all other known carbon phases 
at the given pressure condition fully automatically. We 
can therefore expect that this method can also find with 
high reliability the low energy structures of many other 
materials for which our knowledge is at present still rudi- 
mentary, leading to important advances in the field of 
solid state physics. 
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